Fluctuations of energy flux in a simple dissipative out-of-equilibrium system 
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We report the statistical properties of the fluctuations of the energy flux in an electronic RC circuit 
driven with a stochastic voltage. The fluctuations of the power injected in the circuit are measured 
as a function of the damping rate and the forcing parameters. We show that its distribution exhibits 
a cusp close to zero and two asymmetric exponential tails, the asymmetry being driven by the mean 
dissipation. This simple experiment allows to capture the qualitative features of the energy flux 
distribution observed in more complex dissipative systems. We also show that the large fluctuations 
of injected power averaged on a time lag do not verify the Fluctuation Theorem even for long 
averaging time. This is in contrast with the findings of previous experiments due to their small 
range of explored fluctuation amplitude. The injected power in a system of N components either 
correlated or not is also studied to mimic systems with large number of particles, such as in a dilute 
granular gas. 

PACS numbers: 05.40.-a, 05.70.Ln, 05.10.Gg, 84.30.-r 



I. INTRODUCTION 

Global quantities (i.e., quantities averaged over the 
boundaries or the entire volume of the system) are of 
paramount importance to describe the dynamics and the 
statistics of dissipative out-of-equilibrium systems 
For instance, the injected power is necessary to maintain 
a dissipative system in an out-of-equilibrium regime. A 
system thus reaches a stationary non-equilibrium state 
when a balance between the mean injected power and 
the mean dissipated power is achieved. The usual tools 
of equilibrium statistical mechanics do not apply to such 
systems or are reduced to the Fluctuation-Dissipation 
theorem [2j. Even more, the injected power / is a fluctu- 
ating quantity and cannot be only regarded as a constant 
parameter. Their fluctuations display values that can be 
several times larger than its average, and their statistics 
(even averaged over a macroscopic volume) present large 
deviations [l|, 0, 0] • It is thus of crucial interest to study 
the statistical properties of the energy flux driving a sys- 
tem far from equilibrium, and its relation to its internal 
energy. In some systems, the energy flux fluctuations can 
be directly related to the internal energy by means of the 
Fluctuation Theorem (FT). 

The FT is of fundamental importance for microscopic 
systems far from equilibrium in a stationary state. It 
was first introduced numerically for a fluid under an ex- 
ternal shear Q , then mathematical proof was given @, 0] . 
For a nonequilibrium dissipative system, this theorem 
describes the asymmetry of distribution of a fluctuat- 
ing global quantity I T (energy flux, entropy production 
rate,...) averaged over a time r much larger than its 
typical correlation time r c . For systems close to equi- 
librium or for macroscopic ones, the FT gives a gener- 
alization of the second law of thermodynamics, and also 
implies the Green-Kubo relations for linear transport co- 
efficients when combined with the central limit theorem 



Q. Moreover, it can be applied to nonequilibrium tran- 
sitions between two different equilibrium states leading 
to the so-called Jarzynski equality [9j]. Its derivation re- 
quires the assumption of time reversibility of the system 
dynamics, ergodic consistency, and a certain initial dis- 
tribution of particle states. Finally, it does not require or 
imply that the distribution of time averaged fluctuating 
quantity I T is Gaussian. 

Experimental tests of the fluctuation theorem rela- 
tion have been reported in various systems: in gran- 
ular gases [lOj, in turbulent flows (thermal convection 
[HIE2I] > swirling flows fl3]). in liquid crystals fl4|, with 
an electric dipole [l5[ or a mechanical oscillator |l6j, in 
a two- level atomic system [l^, and by means of a col- 
loid particle [18[ or an RNA molecule [19[ in an optical 
trap. In all these experiments, the fluctuation theorem is 
found to be verified with good accuracy despite some of 
these systems do not satisfy the microscopic reversibility 
hypothesis. Such a good agreement has been also re- 
ported in numerical simulations of granular gases p], 0] , 
turbulence p], [20[ , and earthquakes [l| . The reasons of 
this apparent verification of the FT are two fold: either 
due to the small range of explored fluctuation amplitude 
e = I T /(I) LU, jj, |2(| , or due to the long averaging time 
r needed y, [2l| . Only small relative fluctuation am- 
plitudes (e < 0.8 for r < 20r c ) have been reached in the 
above experiments [TO, EH, E2, EH • Very recently, large 
range of e has been attained, even for r >> r c , by mea- 
suring the fluctuating injected power in an experiment of 
wave turbulence on a fluid surface [22[ • This experiment 
then shows that the FT is not satisfied for high enough 
e. Such a disagreement was also predicted theoretically 
in a system described by a Langevin equation [23[ . Note 
that the breakdown of FT has been recently reported 
numerically [20[ or theoretically [24| in other systems. 

In this paper, the fluctuations of energy flux in an elec- 
tronic circuit are measured to test the fluctuation theo- 
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rem within a large range of accessible value of fluctu- 
ation amplitude (e 3) even for long averaging time 
(t/t c ~ 20). The electronic circuit is a resistor of re- 
sistance R in series with a capacitor of capacitance C 
driven with a stochastic voltage. This circuit can be view 
as an electronic analogue of the Langevin equation [25[ 
which describes usually the brownian motion of a parti- 
cle [26|. It is important to notice that in our experiment 
the dissipation is selected by the system itself. No ad-hoc 
dissipation or thermostat is introduced to ensure the FT 
hypothesis (i.e. the time reversibility of the system). The 
study of the statistical properties of the injected power 
in such a circuit point out three important results: 

First, the probability density function (PDF) of the 
fluctuations of the injected power in the circuit is stud- 
ied as a function of the control parameters (damping 
rate, amplitude of the stochastic forcing). The asymme- 
try is driven by the damping rate: The more the mean 
dissipation increases, the less the negative events of in- 
jected power occur. This electronic circuit is one of the 
simplest system to understand the properties of the en- 
ergy flux fluctuations shared by other dissipative out-of- 
equilibrium systems (such as in granular gases [lflj, wave 
turbulence [22[ and convection [H, [28[ ) . 

Second, we show that the fluctuations of injected power 
averaged on a time r do not verify the fluctuation theo- 
rem at large values of e, even for r >> r c . This occurs 
for values of e larger than the most probable value of 
the injected power PDF. This electronic circuit thus ap- 
pears to be a very useful tool to test fluctuation theorem 
in the different limits of the averaging time and of the 
fluctuation amplitude. 

Third, the injected power in a system constituted by 
an ensemble of N uncorrelated components is then stud- 
ied. This mimics a dissipative multi-component system 
driven out-of-equilibrium without spatial correlation be- 
tween them. The fluctuations of the time averaged in- 
jected power of the N components then verifies the fluc- 
tuation theorem for finite time. This bridges the gap be- 
tween results about the test of the FT for systems with 
low particle number (such as the ones described by the 
Langevin equation), and systems with large number of 
uncorrelated particles (such as in a dilute granular gas). 
The link between them can be understood as a conse- 
quence of the central limit theorem. 

It is well known that electronic circuits are very use- 
ful analogue experiments to study stochastic nonlinear 
problems [29[. However, one could wonder their relevance 
with respect to numerical simulations. Analogue circuits 
get the advantages that any naturally occurring noise 
necessarily has a finite correlation time, and thus avoid 
to pre- select a correlated noise type (Ito-Stratanovivch 
dilemma) in writing the numerical code (29[. Moreover, 
the simulation leads to the accumulation of truncation 
errors, and it takes a longer time to implement and to 
compute [29( | . 

The paper is organized as follows. Section HJ explains 
the experimental setup of the RC circuit. Section [Till 



contains the results about the statistical properties of the 
injected power in the circuit. Some of the experimental 
results of the Section HIJ are then recovered in Section HVl 
with a simple model based on a Langevin equation with a 
Gaussian colored noise (the so-called Orstein-Ulhembeck 
noise) [22|, HtJ • Section [VI contains the experimental test 
of the Fluctuation Theorem for the energy flux in a RC 
circuit. Finally, Section[VI]is devoted to the experimental 
study of the injected power in a system constituted by a 
set of N uncorrelated components, as well as the test of 
the Fluctuation Theorem for its energy flux. 



v(t) 



l(t) ~ C(t)v(t) 
T 



C(t) 



FIG. 1: Scheme of the electronic circuit as an analogue of the 
Langevin equation. 



II. EXPERIMENTAL SETUP 

The experimental setup consists of a resistor, i?, in se- 
ries with a capacitor, C, driven with an external stochas- 
tic voltage £(£) as displayed in Fig. [H The equation of 
continuity for the voltage V across the resistor R leads 
to 



(1) 



where 7 1 = RC. As it will be shown below, the injected 
power in the RC circuit is 



i(t) = 7 v(t)C(t). 



(2) 



The zero mean Gaussian random noise is generated 
by a Spectrum Analyzer (Hewlett-Packard HP 35670 A). 
This noise is low-pass filtered at a cut-off frequency A 
fixed to 5 kHz, unless specified otherwise. The control 
parameter is the noise amplitude D defined by the con- 
stant value of its power spectral density, as an analogy 
to the white noise limit. C is fixed to 1 /iF, and R can 
be varied between 200 ft and 10 kQ leading to values of 
7 between 50 Hz and 10 kHz. The output V(t) of the 
RC circuit is multiplied by the random forcing by 
means of an analog multiplier (Analog Devices AD540). 
The resulting voltage V(t)((t) is proportional to the in- 
jected power (see below) and is acquired with a Digit al- 
to-Analog Acquisition card (AT-MIO-16X) at 100 kHz 
sampling frequency for 10 s, with a precision of 0.3 mV. 
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Equation (pQ) is the analogue of the Langevin equation 
that usually describes the dynamics of a Brownian par- 
ticle of velocity v as [26| 



dv(t) 
dt 



(3) 



where 7 is the inverse of a damping time. / is an ex- 
ternal Gaussian random forcing with zero mean and a 
given autocorrelation function. In the singular limit of 
zero-correlation time (i.e., for a white noise forcing), 
this function reads < f(t)f(t f ) >= foS(t — t') and the 
Fluctuation-Dissipation theorem is satisfied with (y 2 ) = 
fo/(2j), fo being the noise intensity [2[. For a non-zero 
correlation time (as in this study) , the system cannot be 
in equilibrium, and another viscous term different from 
the one of Eq. ([3]) must be used to recover the equilibrium 
state |30|. Multiplying Eq. ^ by v gives 



d 
dt 



v(t) 



21 



f(t)v(t) - ^v{tf 



(4) 



meaning that the energy budget of the system is driven 
by the injected power, f(t)v(t), and the dissipative one, 
^v(t) 2 . This analogy thus shows easily that Eq. ([2]) is the 
injected power in the electronic circuit. 

The aim is now to study the probability distribution 
function (PDF) of the injected power in the RC circuit, 
described by a Langevin equation as the simplest dissi- 
pative system driven out of equilibrium by an external 
force. The objective is to probe the out-of-equilibrium 
statistical properties of the injected power and its rela- 
tion with the fluctuation theorem. It is noteworthy to 
underline that in this simple system the forcing f(t) is 
not in any case a thermal bath. Due to the non-zero cor- 
relation time of the forcing, this system is strongly out 
of equilibrium and the Fluctuation-Dissipation theorem 
does not hold [2]. This is mainly due to the non-Gaussian 
shape of the injected power distribution, in contrast with 
other experimental devices where the injected power fluc- 
tuations are quasi normal [HI, Qi| • 
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FIG. 2: (Color online) Probability density functions of the 
injected power I for two different noise amplitudes: D = 
1.56 x 1(T 3 V r 2 ms /Hz [(a) and (b)], D = 0.75 x 1(T 3 V r 2 ms /Hz 
[(c) and (d)]; and damping rates: 7 = 200 Hz [(a) and (c)], 
7 = 2000 Hz [(b) and (d)]. 
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III. STATISTICAL PROPERTIES OF THE 
INJECTED POWER 

The probability density function of the injected power, 
I(t), is shown in Fig. [2] for different values of the noise 
amplitude D, and the damping rate 7. For all values of 
D and 7, the PDFs exhibit two asymmetric exponential 
tails and a cusp near I ~ 0. Note that this typical PDF 
shape has been also observed in various more complex 
systems (granular gases [Ioj |. wave turbulence [22] and 
convection [12|, |28[). As shown in Fig. [2j the PDF skew- 
ness increases strongly with 7 at a fixed D. Moreover, 
the extremal fluctuations increase strongly with D at a 
fixed 7. 

At a fixed value of 7, the PDFs of I are plotted in Fig. [3] 
for 9 different noise amplitudes. As shown in the inset of 



FIG. 3: (Color online) Probability density functions of in- 
jected power, 7, for D = 0.06 to 1.56 x 10~ 3 V? ms /Hz (see the 
arrow) for 7 = 200 Hz. Inset: Probability density functions 
in the rescaled variable (I — (I))/<ji. 



Fig. 02 all these PDFs collapse on the same curve when 
plotted in the centered-reduced variable, (I — (J))/crj, 
where 07 is the rms value of 7, and (I) its mean value. 
Such a collapse means that all the moments of I scale as 
07. As shown in Fig.[H 07 (as well as (I)) scales linearly 
with D. This linear dependence with D of the PDF of 
I can be recovered by dimensional analysis. Thus, since 
the slopes of the exponential tails scale as when the 
noise amplitude D is doubled, the largest injected power 
fluctuation reached is doubled. 
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FIG. 4: (Color online) Mean (I) and standard deviation ai 
of the injected power as a function of the noise amplitude D. 
7 = 200 Hz. 




FIG. 5: (Color online) Mean (/) and standard deviation ai 
of the injected power as a function of the damping rate 7. 
D = 0.75 x 10" 3 Vrms/Hz. (-): linear best fits of slopes 1.9 
V and 1.59 V, respectively. 



The noise amplitude D is now fixed in order to take 
into account the effect of the damping rate 7 on the in- 
jected power fluctuations. For different values of 7, (I) 
and ai are plotted in Fig. [5l Both moments scale as a 
power law of 7 with two different exponents. Therefore 
no collapse occurs when the PDFs of / are plotted in the 
centered-reduced variable. However, as displayed in Fig. 
[6j both the exponential tails of positive and negative val- 
ues of / show power law dependences with 7. The slope 
of the positive exponential tails scales as ~ ^-i-65±o.05^ 
whereas the negative one scales as ~ ^-i-33±o.05 This 
means that the probability of having negative values of 



injected power decreases faster than the probability of 
having positive ones as the system becomes more and 
more dissipative. 
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FIG. 6: (Color online) Scaling of the PDFs of the negative 
values (left) and the positive values (right) of injected power 
/, for 9 values of D, and 10 values of 7. 

Taking into account both the effect of D and 7, the 
PDF of the positive values of / behaves, far from the 
cusp at / ~ 0, as 

P + (/)~exp(-a+-^y . (5) 

Similarly, the PDF of the negative values of / behaves as 

P_(/)~exp(a_-^X^J ( 6 ) 

where a± are two constants. As shown below in Sect. II V[ 
an explicit formula of the PDF of / can be computed [22| , 
that can capture the properties of the distribution found 
here: a cusp close to zero and asymmetric exponential 
tails (see Sect. [IV)). 

Both D and 7 are now fixed in order to study the effect 
of the random noise cut-off frequency A on (/) and 07. As 
shown in Fig. [7J when A is varied from 3 kHz to 40 kHz, 
the mean injected power is roughly found independent of 
A with our experimental accuracy, whereas 07 scales as 
the square root of A. 

Finally, to summarize all the experimental results, the 
two first moments of injected power behave like 

(/) ~ Dj 1 - 90 and at ~ Dj 1 - 5 ^ - 50 . (7) 

Note that all the previous exponents are experimentally 
measured with an accuracy of ±0.05. Thus, the noise 
amplitude D is found to drive the scale of the injected 
power fluctuations whereas the damping rate 7 controls 
the asymmetry of the PDF of /. 
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FIG. 7: (Color online) Scaling of the mean (I) and standard 
deviation 07 with the cut-off frequency A. (— ): linear best fit 
of slope 0.11 V and 0.56 V, respectively. 



IV. L ANGEVIN-TYPE MODEL WITH AN 
ORSTEIN-ULHEMBECK FORCING 

Using a simple model that has been recently presented 
in Ref. 1221 and will be discussed in details in another 
paper [27[, let us try to recover the above experimental 
results: The shape of the injected power distribution, and 
the scaling of its first cumulants ((/) and 07) with the 
parameters D, 7 and A. 

From Eq. (pp) of the electronic circuit, and the fact that 
the stochastic forcing ((t) is low-pass filtered at frequency 
A, one can write the following coupled linear equations 



dV(t) 
dt 

d((t) 
dt 



jV(t) = 7 C(t) 



(8) 
(9) 



with V(t) the voltage, 7 -1 = RC the damping parame- 
ter, ((t) the colored random forcing, and £(£) the Gaus- 
sian white noise with < >= D5(t — t'), D being 
the noise amplitude. Note that if we only used a Gaus- 
sian white noise in Eq. ([8]), then one find (/) ~ 07 ~ D 
but with no dependence with 7 [2[ contrarily to the ex- 
perimental results [see Eq. (|7J)]. A dependence with 7 
is obtained when using an colored type of noise for 
such as the Orstein-Ulhembeck (OU) one of Eq. ([9]) [2|. 
The colored noise indeed introduces a typical frequency 
needed to simulate the frequency cut-off A experienced 
by the low-passed filtered Gaussian white noise in the 
experiment. 

As shown in Sect. [Ill the injected power in the circuit 
writes I(t) = j((t)V(t). Using the fact that both vari- 
ables V(t) and ((t) are Gaussian with zero mean, the 
PDF(7) can be written in an explicit way [22[. Let us 
rapidly recall the main points of its derivation. First, 



the stationary joint PDF(V,C) writes as a Gaussian bi- 
variate which depends only on the correlation coefficient 
r = ((V) /crycr^ between both random variables [3~0 ], 
where = y/D/{2\) and ay are the rms values of ((t) 
and v(t), respectively. Second, by means of a change of 
variables, the PDF(J = Cy = 1/7) then is computed as 




P(J) = V exp 



7TC 



rl 

c 



(10) 



where c — (1 — r 2 )avcr^, and Kq[-] is the zeroth order 
modified Bessel function of the second kind. One have 
also r = VtTCtTA) [22|, meaning that, at fixed A, r is 
directly related to the damping coefficient 7. Eq. ([TO]) 
then is determined once r is known, i.e. when (/), cry 
and <j£ are known. Since these quantities are measured, 
we can compare the theoretical PDF(J) of Eq. ([TO]) with 
the experimental one with no adjustable parameter. This 
is shown in Fig. [8] for two different values of 7 (or r). The 
computed PDFs display a cusp at / = and exponential 
asymmetrical tails for large values of / in good agree- 
ment with the experimental shapes. As shown in Fig. [HI 
increasing the damping rate 7 leads to PDF more and 
more asymmetrical with less and less negative events. 
The asymmetry then increases when the damping rate 7 
increases. The asymmetry or the skewness of the injected 
power distribution is then controlled by the damping pa- 
rameter 7 (or the correlation coefficient r at fixed cut-off 
frequency A). 




FIG. 8: (Color online) PDFs of // (/): Comparison between 
experiment (— ) and theory [(— • — ) from Eq. ([TU]) ] for two 
different values of the damping rate 7 = 2000 Hz (r = — 

0.45) [black line] and 7 = 200 Hz (r = 0.15) [red (light gray) 
line]. The cut-off frequency A is fixed to 10 kHz. 

For other dissipative out-of-equilibrium systems show- 
ing energy flux fluctuations, an analogue of the parameter 
7 can be found. For instance, in an experiment of wave 
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turbulence on a fluid surface [22[ , the distribution shapes 
of the injected power / by the wavemaker resemble to 
the ones found here: When the fluid used is mercury, the 
PDF(J) is strongly asymmetrical whereas with water, it 
is much more symmetrical. This is due to mean dissipa- 
tion which is different for each fluid. The analogue of the 
7 parameter for wave turbulence experiment is indeed 
related to the inverse of a typical damping time of the 
wavemaker which is linear with the fluid density [22[ • 

With the Langevin-type model of Eqs. ([8]-[9]), one can 
also calculate the first cumulants of I(t). By solving the 
linear part of Eqs. ([8]-[9|), the first cumulants of I(t) in 
the stationary limit read [31] 



a i = 7 



DX 

1+7' 



DA 



AV2 7 l/2- 

In the limit 7/ A << 1, Eq. (JTT]) yields 
(/) ~ D 7 2 A°, 



(11) 



(12) 



(13) 



which does not depend on the cut-off frequency A, and 
Eq. (ED yields 



< D 7 3/2 A 1/2 



(14) 



The range of 7 used experimentally is between 50 and 
2000 Hz, and the frequency cut-off A is in the range from 3 
kHz to 40 kHz. This leads to 7/ A ~ 0.1 in the worst case. 
The first two cumulants of Eqs. ([T3j) and ([T4j) derived 
from the OU process thus are in good agreement with 
the experimental results of Eqs. ©. 



V. RELATION WITH THE FLUCTUATION 
THEOREM 

The smoothing average of the injected power I T is com- 
puted from the previous data of I as 



1 f t+T 

r (t) = - / j(0< 
T Jt 



(15) 



where r stands for the time of average of the signal, which 
is several times the correlation time r c of the injected 
power /. For our experiment, the correlation time r c is 
the inverse of the cut-off frequency, 1/A, which is now 
fixed to 10~ 4 s. 

To describe the asymmetry of time-averaged injected 
power I T distribution, the quantity p(e) is computed as 



lim — In 

t— >co T 



P{e) 



P(-e) 



(16) 



where P(e = I T /(I)) is the probability to have a e equal 
to a certain value I T /(I). p(e) is usually called the asym- 
metrical function [23]. Equation ([T6]) is called the Fluc- 
tuation Theorem which states that, for times r larger 



than r c , this function depends only on e @, 0, Q- In a 
certain limit, Eq. ([T6]) takes the form 



p(e) = pe 



(17) 



where (3 is a dimensionless constant. It means that 
the probability ratio to have a positive value of injected 
power (e) with respect to its negative value (— e) increases 
exponentially with e at large r. Note that a similar re- 
lation called the Gallavotti-Cohen relationship has been 
derived under specific conditions p. The hypotheses for 
deriving Eq. ([T7]) are three: the system should be micro- 
scopically reversible, dissipative and the dynamics on the 
phase space should be chaotic 0, 0, Q • For our dissipa- 
tive system, the reversibility condition is obviously not 
fulfilled. However, let us try to test the relation of Eq. 
([T7|) with our experimental data of injected power. 
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FIG. 9: (Color online) PDF of I T /{I) for various values of 
t/t c = 0, 5, 10, 50, 100 and 200 at a fixed value of 7 = 2000 
Hz. The straight line ( — ) correspond to I T /(I) = and the 
dashed line ( ) to I T = (I). 

Figure [9] displays the PDF of time-averaged injected 
power I T /(I) when r/r c is increased. Several features 
appear. First, the negative injected power events de- 
crease with increasing r until they disappear for r > 5r c . 
Second, when r/r c is increased, the PDF shape for neg- 
ative values of I T /(I) change from an exponential shape 
to a Gaussian one, whereas the exponential shape of the 
positive part is quite robust. Only when r >> r c , the 
PDF shape close to the maximum tends towards a Gaus- 
sian, as one would expect from the central limit theorem. 
In Fig. [9j when r/r c increases, the PDF most probable 
value e* (i.e., where the PDF amplitude is maximum) 
increases slowly from I T /(I) = to 1 (the mean value 
of the injected power). This dependence of e* is shown 
in Fig. [11] as a function of r/r c . This dependence will 
be of fundamental importance when probing the FT (see 
below). 
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FIG. 10: (Color online) Same as Fig. El for 7 = 200 Hz. 



The Large Deviation Function (LDF) /(e) is generally 
denned as 



FIG. 11: (Color online) Most probable value e* of 
PDF(I r /(I}) as a function of r/r c for 7 = 50, 100, 200, 500 
and 1000 Hz. For e < e*, the relation p(e) ~ e is verified at 
finite r, whereas it does not hold for e > e* (see text). 



/( C )= lim -£ln[P( e = I T /(I»], (18) 

T— s>CO 7" 



(19) 



and Eq. ([16]) thus leads to 



p(<0 = /(e) " /("£)• 



The LDF describes the probability of very large and un- 
common events of e. It is consequently very hard to mea- 
sure it. The computed LDF as in Eq. ([T8]) approaches its 
theoretical limit only for large values of r/r c . With our 
experimental data, one can probe large values of r/r c 
and therefore calculate a very accurate estimate of the 
LDF. Developing Eq. (fT9|) up to first order in e, that 
means, regarding only the terms e ~ 0, thus leads easily 
to verify Eqs. (fT6|) and $T7j). This was first conjectured 
by Aumaitre el al. 1] and then predicted in a particular 
system by Farago [23[. But, what would happen if e was 
far from zero? 

The experimental values of the asymmetrical function 
p(e) are shown in Figs. [T2l and [T3l for two different values 
of 7, as a function of e with < e < 3. For small e, 
p(e) increases linearly as expected, then p(e) saturates 
when e is increases further. For each value of r/r c , the 
beginning of the saturation occurs for a critical e value 
called e c < 1. Thus, the linear prediction p(e) ~ e is valid 
at finite r only for e < e c . It is important to notice that 
the saturation value e c of Figs. fT2l and [T3l corresponds to 
the maximum value e* of the PDF (see Fig.fTT]). The fact 
that p(e) 00 e for values of e greater than e c = e* is due 
to the different shapes of the PDF(e) for e < — e* and for 
e > e* (see Figs. l9l and [TO]) . By extension to non finite 
r, this thus means that the FT relation of Eq. ([TT]) does 
not hold for values of energy flux greater than its most 
probable value e*. 



When t/t c increases, it is predicted that the most 
probable value e* of the PDF converges to 1 slowly as 
a power-law of r/r c [24]. This power-law dependence 
is not found experimentally with our data (not shown 
here), and Fig. [TT] clearly shows the convergence of e* to 
1 as a consequence of the convergence of the computed 
/(e) towards the LDF. An analytical prediction for the 
LDF of the injected power distribution has been derived 
for a Langevin equation either with a white noise [23[ or 
with a colored noise (OU) forcing [33]. At high r/r c , the 
shape of our experimental PDF roughly tends towards 
a Gaussian (see Figs. l9l[T0l at r/r c = 200) contrarily to 
the asymmetrical prediction of the LDF with a white or 
a colored noise. However, we have to be careful during 
this comparison due to our low statistics at very long av- 
eraged times (see the vertical range in the Figs. l9l[T0l at 
r/r c = 200). 

Increasing now 7, at fixed r/r c , leads decreasing avail- 
able values of e necessary to probe the Fluctuation The- 
orem (see Figs. fT2l and [T3]) . It comes from the fact that 
when 7 is increased, the number of negative injected 
power events, e < 0, decreases (7 controls the skewness 
of the PDF at a given r c ~ We stress the fact 

that the damping rate 7, and therefore the mean dis- 
sipation, is not chosen in this simple experiment in an 
ad-hoc manner to satisfy time-reversibility. The smooth- 
ing of the signal around (I) also decreased the number of 
available negative events. 

In most of the previous experimental test of the Fluc- 
tuation Theorem 0, El, EI El the limit of Eq. JED 
is well followed, because of the small range of explored 
e < 0.8 at high r/r c < 20. However, very recently, large 
range of e has been measured experimentally [22[ the 
Fluctuation Theorem was not satisfied. In our experi- 
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FIG. 12: (Color online) Asymmetrical function p(e) — 
m [ p(^) ] as a f unc tion of e for different integration times 
t/tc = 1 (★) to 31 (+) at fixed 7 = 100 Hz and D = 1.56 
mV? ms /Hz. 



VI. ON THE INJECTED POWER 
CORRELATIONS 

Dissipative stochastic systems driven out-of-equili- 
brium are generally constituted of several components 
(e.g., in granular gases p], [To[) that may display corre- 
lations in space and time. One can wonder how these 
spatio-temporal correlations change the PDF of the in- 
jected power. Even more, it is important to study their 
relevance in the fulfillment of the Fluctuation Theorem. 

The correlation time of the injected power into our 
simple experimental system can be tuned as a control pa- 
rameter. To wit, the averaged injected power signal I T (t) 
is expressed as a sum of correlated variables where their 
temporal correlations mimics the spatial correlations in 
extended high-dimensional systems (see § VI A|) . One can 
also look at the sum of M independent random variables 
distributed as I(t) (see WIB|) . These two kind of signal 
processing are performed to understand if a set of sta- 
tistically dependent or independent components has an 
effect over the fulfillment of the FT (see JVTC]) . 



A. Correlated components 
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FIG. 13: (Color online) Same as Fig. E] for 7 = 50 Hz. 



For a single electronic circuit, the smoothing average I r 
of the numerically sampled injected power J(t), defined 
in Eq. ([T5]) can be written as the discrete sum over TV 
points, 



1 At) 



1 N 



(20) 



with r = TV At and At the inverse of the sampling fre- 
quency. In our experiment, At is fixed at 10 /is. Since the 
correlation time of the injected power, r c c± 1/A = 100 /is, 
is greater than At, the elements of the sum above have 
a nonzero temporal correlation. 

This smoothing average can be also viewed as a sum 
of N statistically dependent components as 



j T (t) 



N N 



(21) 



k=l 



k=l 



where 1^ (t) corresponds to the injected power of the kth 
correlated component. 



ment, large range of e (up to 3) are also available even 
for high t/t c c± 20. This, thus allows us to test deeply the 
Fluctuation Theorem. As explained above, the FT works 
only for e values smaller than the most probable value e* 
(see Fig. [TT]) . Above this value saturation occurs, due to 
the different behavior of the PDF: for values larger than 
the most probable value, the PDF remains exponential, 
whereas for values smaller than e* it is smoother. Thus, 
large events of injected power are not well described by 
the FT, and lead to the observed saturation of p{e). 



B. Uncorrelated components 

Let us now focus on the case where correlations be- 
tween components are neglected. That is to say each 
component losses its memory of the effect of the rest of 
the system faster than its internal dynamics, such as the 
case of a dilute granular gas where every particle dissi- 
pates its energy by collisions. After each collision, due to 
the low density of the gas, the particle losses its memory 
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of its initial conditions decorrelating the injected power 
events in time. 

We study JV statistically independent variables each 
distributed as I(t). For each time £, each injected power 
Ii(t) of the ith non-correlated component is summed to 
obtain the ensemble average of the injected power, Ij^(t) 
defined as 



(22) 



with Ii(t) distributed as Eq. ([TO]) . This ensemble aver- 
age Ij\r(t) should have different statistical properties than 
the smoothing one I T (t). Indeed, Ij\f(t) results from the 
sum over JV statistically independent components [see 
Eq. ([22]) ]. whereas I T (t) comes from the sum over N sta- 
tistically dependent or correlated components [see Eq. 
GDI. 



C. Results 

The statistical properties of the injected power into 
both systems described above display striking differences. 
Figures [TH show the computed LDFs of I T and Ij\f re- 
spectively defined by ^ In [P(I T /{I))] (see Eq. ([18])) and 
by ±j. In [P(I M /(I))\. These LDFs describe how the fluc- 
tuations of both averages with respect to the mean (I) 
behave when the number of variables taken into account 
in the each sum, JV or r/r c , becomes larger and larger. 
The computed LDF of I T (t) is always asymmetric with 
exponential tails whatever r >> r c , whereas the com- 
puted LDF of Ij\f(t) tends towards a parabola when JV 
increases. 

For the system of statistically dependent components, 
the LDF of the injected power is not parabolic (as it 
should be if its PDF was a Gaussian) as shown in Fig. 
[T4k . The convergence to its asymptotic shape is slow, 
depending strongly on the number of components of the 
system (i.e., of the durations of the time averaging, r/r c ). 
Moreover, when r/r c increases, Fig. [Mh shows also that 
the PDF's maximum slowly tends towards the mean 
value (I), as already noticed in the previous section (Fig. 
[TT]) . As it has been already shown in Fig. [12] in this case 
the Fluctuation Theorem is not satisfied. 

For the JV uncorrelated or statistically independent 
systems, the computed LDF of Ij^ shown in Fig. [Mb has 
exponential tails whatever the value of JV as expected 
from the distribution of Ii(t) [see Eq. ([TO]) ]. When JV is 
increased from 2 to 10, the center of the LDF becomes 
more and more parabolic (its PDF becomes more and 
more Gaussian) as shown in Fig. [T4b . This is due to the 
central limit theorem which can be seen as a quadratic 
expansion of the LDF around the mean of the distribu- 
tion. At higher JV > 10, no deviation from a parabola is 
observed in the inset of Fig. [Mb due to the small fluctu- 
ation values probed. Indeed, if larger fluctuations could 
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FIG. 14: (Color online) (a) Computed large deviation func- 
tions ^ In [P(I T / (I))] of t/t c correlated variables with t/t c = 
3 (□) to 50 (*) for 7 = 100 Hz. (b) Computed large deviation 
functions In [P(Ij^ / (I))] of Af uncorrelated variables with 
Af = 2 (□), 4 (o), 6 (+), 8(*) and 10 (>) for 7 = 100 Hz. 
Inset: same with Af = 10 (>), 20 (o), 30 (+), 40 (*) and 50 
(□). The dashed lines show the mean injected power (J). 



be accessed, one should expect exponential tails in the 
distribution. 

Let us finally test the FT for an ensemble of Af inde- 
pendent variables. The smoothing average of Ij\r(t) over 
a time r is defined as 



1 M rt+T 



(23) 



The asymmetrical function p(Ij\f T / (I)) of the Af indepen- 
dent variables (10 < Af < 100) is plotted in Fig. HH for 10 
different integration times r/r c . Whatever the value of 
JV and t/t c , p(Ij\f T /(I)) increases linearly with Ij^ T /(I). 
Thus, for the system of JV uncorrelated variables (e.g. 



10 




FIG. 15: (Color online) Asymmetrical function p(Itf T /{I)) 
of J\f independent variables, for different integration times 
t/t c = 1 (★) to 31 (+). Af = 10 to 100 with a 10 step. 

a system without spatial or temporal correlations), this 
means that the asymmetric function p(e = Ij\f r /(I)) sat- 
isfies the Fluctuation Theorem as soon as M > 10 [see 
Eq. (fT6|)]. However, one should be careful with this state- 
ment. Our range of accessible fluctuations e decreases 
with increasing J\f and r/r c . Consequently, one can only 
probe the Gaussian part of the PDF(e), but not the ex- 
ponential tails, leading to the linear behavior observed 
for p(e). Larger e values should be reached in order to 
observe the effect of the exponential tails on the validity 
of the FT. 



VII. CONCLUSION 

In conclusion, we have studied the statistical properties 
of the instantaneous injected power I(t) in one of the sim- 
plest dissipative out-of-equilibrium system: an electronic 
RC circuit submitted to a stochastic voltage. The proba- 
bility distribution function (PDF) of I(t) is measured for 
different values of the forcing amplitude and of the damp- 
ing rate 7. It displays a cusp near / ~ and asymmetric 
exponential tails. This typical PDF shape has been ob- 
served in more complex dissipative systems (such as in 
granular gases [10|, wave turbulence [22[ and convection 
jl2l |28|). The relevant parameters of the system can be 
easily changed in our simple experiment. This leads to 
an heuristic understanding of the features of the injected 
power PDF by means of a Langevin-type model. The 
system response V(t) and the forcing ((t) are indeed de- 
scribed by two Orstein-Ulhembeck random variables that 
follow linear coupled Langevin equations [22|, [27[ . Their 



correlation coefficient r = (VQ/crv&C (related directly to 
the mean injected power) is the only control parameter 
driving the asymmetry of the distribution of I(t): The 
larger the damping rate 7, the larger r, and the larger 
the asymmetry of the PDF. Moreover, from this model, 
the scaling of (/) and 07 are found in good agreement 
with the experimental measurements. 

The fluctuation theorem (FT) has then been probed 
by measuring the asymmetrical function p(e) with e = 
I T /(I) : and I T the smoothing average on a time lag r. 
Contrarily to previous experiments, the range of avail- 
able fluctuation amplitude is large (e 3) even for long 
averaging time (r/r c ~ 20). This experiment thus allow 
to probe the FT in the limit of large e and large r/r c . 
We have found out that the FT is only satisfied for val- 
ues of e smaller than the most probable value, e* (i.e. 
the maximum of the PDF of e). For values larger than 
e*, the asymmetrical function is no more linear with e 
but saturates. Thus, the FT does not hold for the large 
available values of e even at large r/r c . This disagreement 
is not particular of this electronic system, but seems to 
be generic to other systems. It has been also recently 
observed with a wave turbulence experiment [22[. This 
model experiment thus appears as a useful tools to probe 
the FT in different limits of averaging time and fluctua- 
tion amplitudes. 

Finally, this electronic experiment can be extend to 
mimic the behavior of a more complex out-of-equilibrium 
systems. To wit, we have studied the injected power 
fluctuations in i) a system of TV statistically independent 
components and ii) a system of M statistically dependent 
components. This latter can be viewed as an archetype 
of a dilute granular gas of uncorrelated particles. The 
Fluctuation Theorem (FT) for the time-averaged injected 
power has then been tested for the case of the correlated 
and uncorrelated systems. In the presence of non-zero 
correlation between components the FT is not satisfied, 
whereas it is satisfied for the uncorrelated system for fi- 
nite average time r. In this last case, the fulfillment of 
the relation is just a consequence of the central limit the- 
orem. Finally, this work also points out that the agree- 
ment with the FT relation is dependent on how the aver- 
aging process is performed (non overlapping bins of du- 
ration t > t c [l0[ or overlapping ones are two different 
processes related to respectively statistically independent 
or dependent components of the system under study). 
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